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I. INTRODUCTION 

The quantum description of the universe is one of the main goals of any quantum gravity approach, since it is 
expected to tame one of the pathologies of General Relativity, the initial cosmological singularity. This is generically 
realized in minisuperspace, i.e. after imposing homogeneity (eventually isotropy) such that the spatial metric tensor 
is described by a fiducial static metric and three (eventually one) time-dependent scale factors. The restriction to 
minisuperspace simplifies the study of the dynamics and it usually allows us to carry on the quantization procedure 
further. 

All these features are realized in Loop Quantum Cosmology (LQC) Bi- LQC adopts some technical tools of Loop 
Quantum Gravity (LQC) in order to quantize phase space variables in minisuperspace. Despite the difficulties 
of LQC to analyze the dynamics of the gravitational field, LQC successes in tracing a singularity-free history for 
the Universe: the Big Bang scenario is replaced by a Bounce occurring when the energy density of the thermal bath 
reaches a critical value. Furthermore, the non-trivial behavior of perturbations in LQC provides potentially testable 
implications in observational cosmology, in particular on the cosmic microwave background radiation spectrum in. 

In |,i we introduced a new model for cosmology in LQC, Quantum Reduced Loop Gravity (QRLG). This model 
is based on implementing in the kinematical Hilbert space of LQC the gauge fixing conditions imposing diagonal 
metric tensor and triads 0- This way we get a reduced Hilbert space, whose elements are based at cuboidal 
graphs and described at links by 17(1) group elements, projected from SU{2) ones. In order to discuss the dynamics, 
some approximations are performed. In particular, the Hamiltonian after the gauge fixing contains some non local 
contributions, which essentially depend on spatial gradients of reduced phase-space variables. We neglect these terms, 
that means we are considering only the “homogeneous” part of the whole dynamics. This approximation is well- 
grounded close to a homogeneous configuration and when the Belinski Lipschitz Kalatnikov conjecture [ll|, [l^ holds. 

Hence, we get a quantum description for the Bianchi I model in terms of the basic structures of LQC, i.e. a 
graph with attached SU(2) quantum numbers, and the dynamics can be analytically investigated. We started such an 
investigation with the easiest case, a three-valent node. We found that the semiclassical limit of the scalar constraint 
acting on such a state reproduces the quantum constraint adopted in LQC [T^ (s ee for earlier attempts 

to relate LQC with LQC via coherent states techniques). Then in [T^ (see alsofl^l we outlined how a consistent 
picture is inferred by describing each homogeneous Bianchi I patches in terms of semiclassical states based at a generic 
graph having N six-valent nodes. In fact, the semiclassical limit of the scalar constraint still reproduces the analogous 
expression in LQC, but the regulator is now related with the number of nodes inside each patch, which is an observable 
on a quantum level (see [l8j|l. 

In this work, we give in details all the technicalities behind the results presented in [l^ . At first, we present a brief 
introduction to LQC in section m and we outline the role of SU(2) recoupling theory in finding the representation 
of holonomy operators. Then, in section imi we discuss the structure of the kinematical Hilbert space in QRLG. We 
define reduced phase-space variables after imposing the diagonal gauges for the metric and the triads and we recall 
how to impose the gauge fixing conditions in the kinematical Hilbert space of LQC, resulting in reduced states based 
at cuboidal graphs with projected group elements at links. We stress how reduced intertwiners arises and how to 
construct the operators corresponding to phase space variables. At this point we present the new technical tool of 
this work, which is the new reduced recoupling rule we adopt to represent the action of holonomy operators. 

In previous papers the recoupling rule for the tensor product of two basis elements (Wigner matrices) 

is defined from SU (2) recoupling theory by projecting twice: on the magnetic numbers of the initial group elements 
and on the spin number of the resulting element. Here, we directly adopt U{1) recoupling theory and we outline how 
it naturally comes out from the full SU{2) recoupling theory when we restrict to reduced group elements. 

In section m we analyze the action of the euclidean scalar constraint operator, which is proportional to the full 
Hamiltonian in the adopted approximation scheme. We start by considering a state based at a three-valent node 
and we find that there are only some minor changes with respect to the result presented in [l^ . However, thanks 
to the new recoupling rule, the extension to states based at six-valent nodes is straightforward and we explicitly 
compute the action of the scalar constraint on such states. In section E semiclassical states are constructed and the 
expectation value of the single-node Hamiltonian is evaluated. At first, we present the case of the three-valent node 
and we point out how we formally get the same result as in [l3l| . but in a more consistent approximation scheme. 
Then, we extend this analysis to a six-valent node. We get the final result for the semiclassical expectation value 
of the full Hamiltonian in section lyll by summing over all the six-valent nodes of the graph. We peak semiclassical 
states around a homogeneous configuration and we define collective variables. The final expression of the semiclassical 
effective Hamiltonian in the large j-limit reproduces the analogous in LQC, with the regulator given by the third 
root of the inverse number of nodes. This is our major result, establishing a fundamental correspondence between 
the effective semiclassical dynamics in QRLG and LQC (and also in Group Field Theory cosmology [iqI - I^ L We 
also discuss quantum correction and we show how inverse-volume corrections are enhanced with respect to LQC, thus 
being a distinctive feature of our model. In section rVHI brief concluding remarks follow. 
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II. LOOP QUANTUM GRAVITY 

The kinematical Hilbert space of QRLG is constructed out of that of LQG by implementing the gauge-fixing 
conditions restricting to a diagonal metric tensor and diagonal triads. In LQG the kinematical Hilbert space T-L is 
constructed by quantizing gravity phase-space parametrized by holonomies and fluxes and it finally reads 

n = , ( 1 ) 

being the gauge-invariant space based at each graph T, i.e. it is the space of gauge-invariant functions of L-copies 
of the SU(2) group, L being the total number of links within T. The direct sum in ([1]) extends over all the piece-wise 
analytic graphs in the spatial manifold. Basis elements in T-L are invariant spin networks, whose explicit expression 
reads 


>= n n ( 2 ) 

rier i 

in which U-l' [hi) are SU(2) irreducible represetations (Wigner matrices) associated with holonomies hi along each link 
I G r, while Xn denote the invariant intertwiners implementing SU(2) invariance at each node n. The products in Q 
extend over the nodes n G T and all the links I emanating from each n. The symbol • means the proper contraction 
between the indexes of intertwiners and Wigner matrices. 

The action of the relevant operators is defined according with canonical quantization. The fluxes across the dual 
surfaces to a link I act as the left- or right-invariant vector fields of the SU(2) group based at L As soon as functions 
of holonomies are concerned, their action is that of multiplicative operators. The representation in the spin network 
basis is given by expanding first the considered operator in SU(2) irreducible representations and, then, computing 
the product of each term with the Wigner matrices of the basis elements. This procedure can be performed once a 
prescription for the product of SU(2) irreps is given; this is SU(2) recoupling theory. The basic law of SU(2) recoupling 
theory can be represented graphically as follows 



( 3 ) 


where the triangles denote SU(2) group elements, the lines are Kronecker 5’s over magnetic indexes and the three- 
valent nodes represent Clebsch-Gordan coefficients. The sum extends over all the admissible representations fc, i.e. 
\jl -j2\ <j< jl +j2- 


III. KINEMATICAL HILBERT SPACE IN QRLG 

In QRLG one chooses the reference frame in which the spatial metric tensor is diagonal along some fiducial directions 
parametrized by the coordinates x* {i = I, 2, 3): 

ds^=N^dt^-df df = al{dx^f + al{dx^f + al{dx^ f , (4) 

M being the lapse function while ai denotes the scale factors, and one fixes inverse densitized triads such that they 
are diagonal too. 


= bi 1=^0 “203 \p2\=(la3,ai b3|=^o«ia2, (5) 

Vb = ^0 being the fiducial volume of chosen space region. Other phase space variables are the diagonal components 
of connections Cj = which in terms of the scale factors read 


Ci = 



( 6 ) 


where ' denotes time derivative. 

The gauge-fixing conditions are realized in the kinematical Hilbert space of LQG by restricting the class of admissible 
graphs to cuboidal ones, having links along the fiducial directions, and by a proper projection of SU(2) group elements. 
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This projection is realized such that at each link k along a;* it is attached a representation of the U(l) group obtained 
by stabilizing the SU(2) group along the internal direction ui = Ui, where 

Ml = (1,0,0) M 2 = (0,1,0) M 3 = (0,0,1). (7) 

This is implemented by projecting the magnetic indexes of the Wigner matrices inside the expression ([2|) for SU(2) 
spin networks on the SU(2) coherent states | ± ji,ui) having maximum or minimum magnetic number along ui. 
Graphically, we can represent the new basis elements by replacing the Wigner matrices at each link I as follows 


(j, n\m, ui)(m, ui\D^ ih)\m, ui){m, ui\j, r) = 


3 






m = ±j, ( 8 ) 


in which the matrix Ri denotes the rotation mapping the direction m; into M 3 and the breaks means the projections 
on the maximum or minimum magnetic numbers. Hence, basis elements in QRLG on a given cuboidal graph T can 
be written as 

(h|r,mi,x„) = ]^(ji,x„|mi,u;)]^ mi = ±ji (9) 

n^r I 

where the products Oner 0; extend over all the nodes n € T and over all the links I emanating from n, 
respectively, while there is no contraction of indexes since all the quantities appearing in ([S]) are complex numbers. It 
is worth noting how some nontrivial coefficients are induced at nodes, (ji, x„|mi, u/), which behaves as one-dimensional 
intertwiners in QRLG. The space spanned by (0 is denoted by 

Flux operators can be inferred by projecting those of LQG and it turns out that they are nonvanishing only when 
they are smeared along the associated dual fiducial surfaces, i.e. ^Ei{S^) = 0 for i ^ j, and they just read the 
magnetic index of each projected representation when the the surface and the link are dual: 

n 0- (10) 

The investigation on the action of holonomy operators requires to fix a proper recoupling theory. In we choose 
the following recoupling rule 



in which ni = ±ji,n 2 = ±^2 are the magnetic quantum numbers and a second projection on the reduced Hilbert 
space forces, in the coupling channel k, the magnetic number ni -|-n 2 to be maximum or minimum, i.e. k = |ni -I-M 2 I. 
In this work, we replace the recoupling theory above with the following one 



which means that we are now recoupling only projected Wigner matrices (and not lines). The two recoupling rules 
dm and (| 12 l) differ when we consider two group elements, based at the same link, which are projected one into the 
maximum and the other into the minimum magnetic components. The new recoupling theory can be obtained directly 
from ©: when we replace SU(2) with 17(1) group elements into triangles, the group element on the right-hand side 
does not depend on the spin k and the sum just provides a resolution of the identity leading to (HU). The formula 
(HU implies then that the reduction commutes with the product operation on the spinnetwork states, namely if we 
denote as tpi{gi) tp 2 {gi) two gauge invariant spinnetworks (having their own intertwiners but the same group element 
associated with the link 1) and with P the projection operator, we have P[ipi{gi)] P[ip 2 {gi)] = P['tpi{gi)i^ 2 {gi)]- Note 
that this property does not hold with (ITTI) . 

In the following we will outline how the new choice gives us an expression of the scalar constraint operator matrix 
elements, which simplifies the analysis of the semiclassical limit, allowing us to consider the generic case of six-valent 
nodes. 
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IV. ACTION OF THE SCALAR CONSTRAINT OPERATOR 


As discussed in [I^ , we consider only that part of the scalar constraint whose Poisson action preserves the diagonal 
gauge condition for the metric. This means that we are neglecting the following contributions: i) the scalar curvature 
of the spatial manifold, ii) non-local terms which arise because of the gauge-fixing conditions. Both these terms can 
be neglected close to homogeneity and in the limit in which the BKL condition holds for the inhomogeneous extension 
of Bianchi I model. 

Within this scheme, the scalar constraint is proportional to the Euclidean part of the full scalar constraint He, 
and we quantize it by taking the expression of LQG [2^ - 1^ (regularized via a cubulation of the spatial manifold) and 
replacing the operators in TL with those in ^TL. 

Hence, we study the following Hamiltonian 


= 


1 


IbTrGy' 


■J2^h%[n] = 


1 




EE 


^ \N], 


(13) 


where the summations extend over all the nodes n of the graph the state is based at and over all the triples of links 
emanating from the same node, while the operator [N] acting on n-valent nodes is given by 

E O 


R Tjn 


[TV] := 


Ai 


SiT^lp c{n) 


Ar(n) Tr c(n) = 2 


-)n—3 


(14) 


A/'(n) being the lapse function in n, while the holonomies are in the fundamental representation. The trace now 
means summing over the two possible projections into we can perform: on the maximum {fi — ^) and on the 
minimum (p. = — i) magnetic numbers. Note the use of a different normalization factor with respect to (22l - [^ . 
due to a regularization based on a cubulation instead of a triangulation (there is a factor 6 between the volumes of 
a parallelepiped and of a tetrahedron sharing three converging edges). The coefficient c(n) is the total number of 
non-coplanar triples. 

It is worth noting how the volume operator is diagonal, such that the matrix elements of ^TL in the basis (jH]) 
can be analytically determined. 

Let us first consider the three-valent node discussed in 



where we called the three fiducial direction x, y, z and we put a loop in the plane xy. For simplicity, we assumed that 
the Wigner matrices in the state |n^)/j have all been projected on the maximum magnetic number, thus to/ = j/. 

The operator ()13p will generically attach loops in all the planes, but since we are considering a non-graph changing 
operator we are going to discuss only the case in which it provides a loop overlapping the one already present in |n^)fl. 
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The action of ^TL coincides with the one discussed in [l3l| , until we have to recouple reduced basis elements and use 
the new recoupling rule (1121) . Hence, looking only at the central node we have 



If one now recouples the U(l) group elements hx and hy one ends up with the following expression 


Tr 


R -L—lRCrR 








j. 1 




It worth noting how the state above does not look like an elements of since the coefficient at the central node is 
not an intertwiner. Indeed, since intertwiners are just numbers, the state above is proportional to a basis element in 
^T-L with a factor given by the ratio of the coefficient at the central node and the right intertwiner. One must just 
check the right intertwiner does not vanish and it is the case here provided that the original intertwiner in (jlSp is 
nonvanishing (one must just verify that the triangular inequality holds for the involved spin numbers). 

By repeating the same analysis for the other (planar) nodes, one ends up with the following expression: 
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Tr 




h:rV^hs^\\n^)R = 

= Y \/jx jy {jz + m) s(/^)C'i i 1-1 

M=±3 

Sz I 



(18) 


The coefficient left from the intertwiner reduction at the node are just standard SU (2) intertwiners contracted with 
fixed rotation matrices. They can be recoupled node by node with the original intertwiners, producing a 9j symbol 
at the central node and 6j symbols at the planar nodes summed over the whole tower of intermediate representations 
allowed by the tensor product of the intertwiners. This sum is the main difference with the results of [l3| where only 
the term arising from the tensor product with maximum total spin was taken into account. 

However, it is convenient to keep the products of single intertwiners as in the previous formula, since it outlines 
how the action of the Hamiltonian does not depend on the valence of the node, thus it can be easily extended to 
more-than-three-valent nodes. Given a six-valent node in a cuboidal lattice 



(19) 
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where the position of the node is described by a vector n, with i = x,y,z denotes the vectors along the fiducial 
directions connecting with the nearest nodes, and the links jft^i with i = x,y, z are outgoing from the node in n, we 
get the analogous of the expression (IT^ 


Tr 


/is. ^ /is.J |no,o.o)fl = 



( 20 ) 

where we consider just one node in the position n = (0, 0,0). This result has then to be summed over all the possible 
permutations of the links within the chosen triple. 


V. SEMICLASSICAL LIMIT OF THE HAMILTONIAN 

The definition of coherent states in QRLG can be performed accor ding with the tools adopted in full LQG 
In particular, we define semiclassical states based at a graph T as in |1^ : 

m, ner ier 

where the functions are determined by iJ/ = ‘ belonging to the complexification of the SU(2) 

group and giving the classical values for the holonomy hi and the dual flux Ei around which the state is peaked. The 
definition of these classical quantities involves the physical varialbes ci and pi for £q = e, e being the coordinate length 
of the link /, while the coordinate area of the dual surface is e^. The explicit expression of ^/’^/(to;) reads 

= + ji = \mi\. (22) 

The insertion of such functions insures good peakedness properties in the limit ^ 1, since they approach 

gaussian functions, i.e. 

^fi,imi) ^ i2ji + , 


(23) 
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times some phase terms which are representations of the classical holonomy For a constant 

physical connection [? ] c/ along the link I, the parameter 0i equals ±c/ and the sign depends on the relative orientation 
between the link and fiducial directions. Similarly, the physical momenta equals the flux = STT^lpfhi, rhi being 
the center of the gaussian function in (I23p . We can now construct coherent states associated with the three-valent 
dressed node (HID, which reads as follows 
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Thanks to (flSl) the expectation value of the scalar constraint on such state is given by 
g n^) = -4*(87r7;|)^Ar(n) ^ 



\ijx jy {jz +Ai) 1 _ 1 

3^,3v,3z,jl ,/Ja,,Ma=±5 




My 






J?. 




3z 1 

Uz) 


r 

'^Hi {jx + Mx) Oy ~ ^^y) 


jx +. 

H. 






JZ+My 

(j/ + My) '^H, (jl + Mx) 

<"y ^ lx 

Ji+My Jz+Mx 


V 


r'‘XX+'‘« 


1 

'^H^ijz) 



3x /\3y 
R; 




'^Huijx) ^H.ijy) 






/ V 


/ 


(25) 


The expression above can be simplified using the fact that the intertwiners appearing should be all normalized. 
This can be done by dividing each three valent node by 



In the coherent states (1241) . this normalization must be done twice: for both the intertwiners in the basis elements 
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and the intertwiners in the coefficients (since the latter are dual to the former). This corresponds to use a normalized 
intertwiner basis, for which each intertwiner is just a phase and the expression above is equal to 1. Having normalized 
intertwiners, the full state n) is normalized too, be. 

=1. (27) 


and (ED) becomes 




(j.)«'ffz, (j.) 


0 



Ux+ 


E E E E 

/-1-1 ,, ,,-Ul -Ul 


jy Uz + m) 


_1 

2 




Hi,. 


I-^X 



{jl + (jl) 


'^HiyUv f^y)'^Hiy{jy) 






This expression coincides with the leading order term of Eq.(71) in [l3l| modulo the coefficient, but here no large j’s 
limit has been performed yet. 

Hence, we can use our previous results to evaluate (05)) . In particular, taking the variances ai = with fc > 1 

[^ . we can perform a saddle point expansion around the expectation values ji of the functions '^/h (which are almost 
gaussian), while the corrections are order , thus negligible in the large j’s limit. 

This way, one gets 


{'i'H n ^)^ jx jy {jz + A*) s(m) sin sin, (29) 

y.=±l/2 

and the only differences with Eq.(74) in [l3l| are that we are now using physical phase variables and we have a different 
coefficient. We can expand the square root \/jz + for jz ^ 1/2, so getting 


VJT+li s{fi) 


1 

jz 


1 + 


1 


(8j.)^ 


(30) 


fi=±l/2 

and by inserting physical momenta pi for a single cell (smeared on the dual surfaces to the links U) one gets 


{'i’H \'Sh n^) ~ sincrcsincy . (31) 

It is worth noting how such a new large j’s expansion provides the so-called inverse-volume corrections, which are 
of 1/j order, thus generically bigger than those related with the saddle point approximation for gaussians. 

The extension to the six-valent case is straightforward. The action of the scalar constraint is the same as in the 
three-valent case. The only difference is that now one has to use the normalization of six-valent reduced intertwiners: 
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obtaining 


no.o.ol^i?" g no.o,o) =^ XI XI X 

Ma:,My==l=^ =l=§ 3x-,jyJz,jSx>y'0e 


(32) 


(jz) 





Ux + (jx) 4/^ {jy — IJ.y)'i’Hi^ (jy) 


jx jy (jz + yj) s{p)Ci 1 1 _i 
2 2 2 2 


-1 


Oex,y + {je:c,y) ^ Hi^ {jey,x ^^x)^ Hi {jey,x) 

‘■ex,y ^ ^x,y ey,x ey,x 

X 


(33) 


This expression coincides with (l28l) and the final result of the single-node Hamiltonian is given by (EH summed over 
ordered permutations of links within a given triple and over all triples around n. On a local homogeneous configuration 
(j.e. when the expectation values of spins along a given direction does not change across the node) the first sum just 
provides twice the sum over ordered permutations x ^ y ^ z, while the second simply cancel the factor c(n) in d. 
Finally, the expectation value of the single-node Hamiltonian is given by 


i'i’H non,o| X *^ 0 , 0 , 0 ) 




sin Cx sin Cy 


Px Py 


sin Cx sin Cy 


Px Py 


Pz 


sin Cx sin Cy 


(34) 


VI. COLLECTIVE VARIABLES 

Let us consider each Bianchi I local patch filled with N = N^NyNz six-valent nodes, Ni being the number of nodes 
on a dual fiducial plane. We want to define the phase-space variables {Ci, Pj} describing each local patch as a whole. 
The semiclassical state is constructed by imposing homogeneity for semiclassical expectation values, such {ci,pj} are 
the same for all nodes and the collective (semiclassical) variables are obtained via the following rescaling 

Px = NyNzPx Py = NzNxPy Pz = NxNyPz , (35) 

Cx = NxC, Cy = NyCy C, = V, c,. (36) 

It is worth noting how such a rescaling is due to the fact that we are using physical variables, whose definitions contain 
the length or the area of the fiducial cell. 
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The expectation value of the scalar constraint is the sum over all nodes of the single-node one (1341) and, in terms 
of collective variables, it reads 


2 




N,N, 


X ‘■'y ‘ 


sin^sin^ + NyN,^ 


IPyPz . Cy . C, 

— sin — sin-h Nz 

Px Ny N, 


IPzPx . C, . C, 

—=— sm-sm — 

Py N, 


(37) 


This expression formally coincide with the quantum operator adopted in LQC for the Bianchi I model [2y, in 
which one takes as regulators the inverse of the number of nodes on fiducial planes, i.e. 




(38) 


Therefore, we infer a correspondence between the quantum dynamics of LQC and the semiclassical description of 
QRLG, similarly to the findings of Group Field Theory cosmology [2^[2l|. 

The regulator is due to the presence of several nodes inside each local patch. The choice of the regularization 
schemes in LQC correspond to the choice of the regulator as a function of the momenta variables. This means that 
the regularization scheme is fixed once we fix the dependence of the number of nodes on Pi. 

At the same time, the number of nodes together with spin quantum numbers enters the definition of the collective 
momenta variables Pi as follows 


(39) 

Henceforth, the ambiguity in the regularization is due to the fact that we cannot uniquely fix Ni as a function of Pi, 
because another fundamental quantity (ji) is involved. For instance, the old regularization scheme is reproduced for 
a constant number of nodes, such that Pi is proportional to ji, while the improved scheme can be inferred (with a 
bigger A) for constant spin numbers [l^ . 

The study of corrections reveals that we get different predictions with respect to LQC; while holonomy corrections 
stay the same, inverse-volume ones differ. This can be seen from the expansion pop . which provides the following 
corrections to the term in (l3T]) 


1 

Vfi 


1 

VFi 



which, in turn, corresponds to the following modification in (1371) 


1 1 / Tr^lp\ 


(40) 


(41) 


The factor 
corrections 
compatible 


N multiplying the next-to-leading order term implies that we get an enhancement of inverse-volume 
with respect to LQC. As consequence, the set of admissible regulators is restricted with respect to that 
with the lattice refinement scheme of LQC 331 (see [l^ for details). 


VII. CONCLUSIONS 

We analyzed the action of the euclidean scalar constraint operator on six-valent nodes in QRLG. This allowed us to 
investigate a realistic model for the quantum universe described in terms of a discrete graph, whose building blocks are 
cuboidal cells. We changed the recoupling rule with respect to previous works @,[13 and we constructed semiclassical 
states, such that the expectation value of the Hamiltonian equals the classical Hamiltonian in the large j limit. This 
is a consistency check on the viability of our model and we recognize the new recoupling rule (I12p as the correct one. 

We also derived the effective Hamiltonian of LQC (1571) . in which the regulator is given by the total number of nodes 
along fiducial directions within each homogeneous patch. The same result has been obtained in the framework of 
lattice refinement and in GFT cosmology [l^, [2l| . We briefly discussed how different regularization prescriptions 
can be inferred in our model by assuming that the number of nodes is a function of semiclassical variables. However, 
the actual implementation of this idea would imply a new definition of the scalar constraint operator, which is behind 
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the scope of the present work. Our analysis is restricted to the non graph-changing version of Thiemann’s Hamiltonian, 
which can be seen as the most conservative definition for the quantum dynamics. Other proposal can be considered, 
like for instance the one in [s^. In this respect, QRLG is a privileged arena to test the phenomenological implications 
of different definitions of the scalar constraint on a quantum level. 

A distinctive feature of our model with respect to LQC is the magnitude of inverse-volume corrections. These 
corrections are enhanced, thus they could have a stronger impact on the cosmological observables. In particular, they 
could significantly affect the dynamics of perturbations and potentially induce sensible modifications to the cosmic 
microwave background radiation spectrum. For this reason, the study of perturbations is the most urgent issue of 
QRLG. However, it requires a great theoretical effort, since the dynamics of perturbations is described by those term 
in the Hamiltonian which are neglected in the homogeneous (or BKL) limit. A first step in this direction is to consider 
a non-vanishing scalar curvature, adopting the quantization procedure described in [sH, |3^. 

Future developments will also deal with the introduction of matter fields in the quantum universe described by a 
collection of six-valent nodes. In this respect, the tools of loop quantization are to be adapted to QRLG framework. 
If successful, this analysis will provide us with both the effective dynamics of field on a fixed quantum background 
and the effective contribution to the right-hand side of Einstein equation (within the adopted approximation scheme). 
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